Experimental investigation of the freely cooling granular gas 
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Using diamagnetically levitated particles we investigate the dynamics of the freely cooling granular 
gas. At early times we find good agreement with Haff's law, where the time scale for particle 
collisions can be determined from independent measurements. At late times, clustering of particles 
occurs. This can be included in a Haff-like description taking into account the decreasing number 
of free particles. With this a good description of the data is possible over the whole time range. 



In the study of granular media it has proven useful to 
define states on the basis of the solid, liquid and gaseous 
states of a molecular substance [lj. In this vein of defini- 
tion, a granular gas is a dilute system of macroscopic par- 
ticles in random, quasi-Brownian motion. In contrast to 
a molecular gas, collisions between granular particles are 
only partially elastic. Hence, energy has to be constantly 
injected into the system and one expects the manner of 
excitation to have a non-negligible effect on the state of 
the system Q. However, the independent characteristics 
of the granular gas should be governed by the statistical 
properties of collisions and by energy losses due to inelas- 
ticity. As the behaviour of an excited granular gas can be 
divided into short-time incipient cooling states between 
single excitation events, an investigation of the cooling 
process allows us to study fundamental characteristics of 
the granular gas independent from the specific manner of 
excitation @,H, 

Haff [fl derived a hydrodynamic theory of granular mo- 
tion. The resulting cooling law states that the kinetic en- 
ergy of a spatially isotropic granular gas without external 
driving should decrease like 1/t. This behaviour emerges 
after a characteristic time determined by the density, ini- 
tial speed and particle properties. 

An experimental study of the cooling of a 2-d granular 
gas on a surface together with a comparative MD simula- 
tion has been reported in [7| , with a special emphasis on 
clustering behaviour. Haff's law was not observed, which 
might have been due to additional energy loss by surface 
friction. 

To extend Haff's picture, the freely cooling granular 
gas has been investigated in several analytical and sim- 
ulation studies. Here the main issue is to describe the 
process in the presence of inhomogeneities, presumably 
introduced by inelastic collapse. The authors incorpo- 
rated clustering and inelastic collapse [1, [j| a s well as 
a velocity-dependent restitution coefficient [TO, QjJ . Al- 
though all studies agree in a decrease of cooling with 
time, the resulting cooling exponents are still a matter of 
debate. In all studies the quantity studied is the granu- 
lar temperature (v 2 ), which is, however, increasingly ill- 
defined for clustering states. Furthermore, as we will see 
below in a simple derivation of Haff's law, the speed of the 
particles is a more fundamental quantity. Also, the cool- 
ing process mainly involves the non-clustered particles, 
and in a video microscopy setup with a limited spatial 
resolution it is practically impossible to determine any 



very small particle speeds in the clustered parts of the 
sample. Therefore, the observed quantity in our study is 
the mean particle speed outside of a possible cluster. 




Figure 1: Cooling process for mechanically (left) and magnet- 
ically (right) excited system. Snapshots taken at the driving 
switch-off (a), the Haff time (b) and after cooling has com- 
pleted (c). 

We will briefly discuss the physics behind Haff's law. 
Assume a granular gas in motion, but with no exter- 
nal forcing. If the density and momentum are spatially 
isotropic, the time evolution of the energy T, which is 
purely kinetic, can be described by the loss due to inelas- 
tic collisions taking place at the rate r _1 : 
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Here we define the restitution coefficient p as the amount 
of kinetic energy retained after a collision. The collision 
time r is given by s/v, where v is the speed and s the 
mean free path. In the dilute limit, s is given by the 
density no and the cross section of the particles cr, s = 
l/(no<r). As T ex v 2 , equation (pQ) can be transformed to 
one describing v, the integration of which yields 
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which is Haff's law in its original form, where the Haff 
time th specifies the system relaxation time scale and 
depends only on the restitution coefficient g and the par- 
ticle size. 

As will be discussed below, the assumption of unifor- 
mity does not hold for a system cooled down to a cluster- 
ing state, which is of course spatially correlated. Thus the 
mean free path of particles outside the cluster increases 
and we expect the cooling exponent to deviate from — 1 
for large times and highly clustered states [l2|, Qjjl, [lil ]. 
For this, we assume a collision time r = l/(n(t)av) de- 
pending on the time dependent number density of par- 
ticles outside the cluster, n(t) = nof(t). Inserting this 
into (pQ) and integrating leads to: 
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Note that for constant particle numbers this is equivalent 
to Haff's original law. f(t) can be determined experi- 
mentally (see Fig. [2j) and integrated numerically from 
the measured data, such that a parameterless fit to the 
cooling data is possible. 

The above description will hold until the mean free 
path has grown to the size of the container. At this 
point, which roughly corresponds to the state where only 
few free particles are left, we will return to a Haff-like 
behaviour, where the container size has to be used as the 
mean free path and the restitution coefficient of collisions 
with the walls applies. 

In order to observe the cooling process experimentally 
in three dimensions, we have to compensate gravity, oth- 
erwise the collision statistics would be dominated by col- 
lisions with the bottom of the sample cell after an excita- 
tion switch-off [Hj]. While this compensation is p ossible 
in principle with parabola flights or satellites [16[ , it can 
be accessed more easily by the use of diamagnetic levita- 
tion, which is highly controllable and reproducible. Dia- 
magnetic levitation has previously been applied success- 
fully in studies of granular demixing [l7| , but we are not 
aware of applications in the field of fundamental proper- 
ties like velocity distributions and cooling behaviour. 

A diamagnetic particle in an external magnetic field 
B acquires an induced magnetic moment antiparallel 
and proportional to B and thereby a potential energy 
U oc B 2 . Accordingly, in an inhomogeneous field, the 
force acting on the material is determined by the prod- 
uct of field and field gradient, F oc B V£?, as can be 




time (s) 

Figure 2: Number of particles outside the cluster vs. time 



easily verified by taking the spatial derivative of the mag- 
netic energy. If the resulting ^-acceleration compensates 
gravity, 

a z = ^B-d z B = g, 

the particle can be held in a stable state of suspension 
comparable to actual weightlessness [l8|, [l9| . 




CCD camera 

Figure 3: Setup for magnetic shaking, including sample cell, 
field modification coil, diffuse illumination and CCD video 
observation. 

Our setup consists of a helium cooled superconduct- 
ing magnet coil (NbsSn, Oxford) with a cylindrical room 
temperature bore of 4 cm diameter which is accessible 
for experiments. The coil can generate field strengths up 
to 20 T, which corresponds to a maximum B z -d z B z of 
approx. 1800 T 2 /m. Fig. [H shows the potential shapes 
calculated from the magnetic field chart supplied by the 
coil manufacturers. As can be seen, the residual accel- 
eration in a volume around 1 cm at the levitation point 
(r, z = 0) are of the order of a few 10 -3 g, so that we can 
effectively speak of milligravity. The radial potential is 
bowl-shaped, the height can be tuned by changing the 
maximum coil field Bq for a given type of diamagnetic 
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material. It is not possible to flatten the radial potential 
entirely. 
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Figure 4: Potential maps (gravitational potential subtracted) 
in z and r for the magnet used in our experiment. 

The levitation point is determined by field strength 
Bo and geometry. Hence it is possible to shift this point 
periodically by modulating the field by an amount AB 
using an additional coil, effectively shaking the particles 
magnetically. The resulting acceleration is proportional 
to BqAB and acts on all particles in the sample in a 
similar manner, with variations due to particle height 
and field inhomogeneities. This spatially almost uniform 
excitation provides a certain amount of chaos and even- 
tually leads to a decrease in collective sample movement 
in favour of random single particle motion. Any remain- 
ing collective motion can be compensated by choosing a 
centre-of-gravity reference frame. 

We also use the common method of exciting a gran- 
ular gas, which is to put the sample on top of a me- 
chanical shaker, e. g. a loudspeaker. In the case of a 
levitated sample, the moving piston barely touches the 
lower end of the sample. Particles moving downwards 
to the topmost piston position are re-injected upwards 
with a comparatively high energy, which is dissipated by 
chaotic collisions with the bulk particles. It is useful to 
have access to two fundamentally different states of ex- 
citation: we are able to examine how far the excitation 
method characterises the continuously driven gas and to 
what extent different initial states influence the cooling 
process. 

To get rid of air friction and hydrodynamic interactions 
as well as to avoid water condensation on the particle 
surfaces, the sample cells are built to hold a weak vacuum 
for several hours after depressurisation. 

A schematic of the setup for the magnetically shaken 
case is depicted in Fig. 03 The sample cell is positioned 
at the point of maximum field gradient located near the 
top end of the superconducting coil. The small copper 
modification coil is driven with a square wave current 
from a computer controlled AC power supply. The setup 



is illuminated from above by a diode array and a Teflon 
diffuse layer and observed with a CCD camera from be- 
low. In the case of the mechanical excitation, a plastic 
rod connects the sample cell downwards to a commercial 
speaker membrane driven by the power supply mentioned 
above. In this case the setup is illuminated from below 
and observed from above. 

The system's ^-projection is observed with a lipstick 
camera and further processed frame by frame by a blob 
recognition software (RSI IDL). Thus, we record velocity 
snapshots by measuring the distance between next neigh- 
bours on consecutive frames at a constant frame rate of 
120 fps. Although there is a risk of misidentification and 
thus a systematic underestimation of the velocities of fast 
particles, this method is reasonably correct for the rela- 
tively small speeds relevant in cooling processes. 

In the magnetically shaken case, our sample consists 
of approx. 50 round Bismuth shots weighing 2 ± 0.1 mg 
each; in the speaker setup the particle number is ap- 
prox. 90. Bismuth, which is highly diamagnetic, requires 
a field of 13.7 T (with B z ■ d z B « 1000 T 2 /m) to levitate. 
The metallic nature of the sample prevents electrostatic 
charges, so there is no additional Coulomb interaction. 
However, Bismuth has an extremely low restitution coef- 
ficient of approx. 0.1-0.2, so we are operating in the limit 
of a highly dissipative granular gas. 

The energy input from the copper coil is very small, so 
that the magnetically excited system will only produce 
reasonable particle speeds if the shaking is done at a reso- 
nance frequency determined by the system geometry and 
the potential shape. Thus, both systems are excited at 
a non- varying frequency of 1.57 Hz. In a measurement 
series, the system reaches an excited steady state over a 
heating period of 10 s in the mechanically and 20 s in the 
magnetically shaken case, after which the speaker or coil 
is switched off. Video recording starts precisely at the 
onset of shaking and is terminated manually when the 
system has relaxed into a clustered state. Each experi- 
mental series consists of 50 such movies, the well-defined 
starting time of the movie allows for an easy determina- 
tion of the onset of cooling to in each case, so we do not 
expect significant systematic errors in averaging over the 
whole series. 

In Fig. [5] we have plotted the mean speed for the re- 
maining free particles outside the cluster. The experi- 
mental data averaged over 50 experiments are shown in 
light grey, the darker line includes an additional smooth- 
ing in time over 50 frames each. In all plots, dotted 
lines delineate Haff's law (Eq. [2|) and dashed lines the 
corrections for time-dependent particle numbers (Eq. [3]). 
All parameters are known from independent experiments: 
the excited state yields vq, while the container size, par- 
ticle number no and size a are fixed. Thus there are 
no adjustable parameters in the curves. The cooling be- 
haviour starts to deviate from Haff's law as soon as the 
cluster begins to form. This can happen at later times 
than the typical Haff time, as can be seen in the me- 
chanically shaken sample. In the magnetically shaken 
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Figure 5: Mean particle velocity vs. time for the mechanically (left) and magnetically (right) excited system. Haff's law 
without (solid, see Eq. [2} and with (dashed, see Eq. [3} cluster correction. Time starts from end of shaking, arrows mark the 
Haff time th 



sample, this roughly coincides with the Haff time, which 
is marked by an arrow in Fig. [5l See also the snapshots 
in Fig. [H The deviations are in very good agreement 
with our calculations derived from the extracted particle 
numbers. Note also that we have no free fitting param- 
eters left, as speed and mean free path can be extracted 
from the video data and the restitution coefficient can be 
established independently. 

Especially in the mechanically excited system, we can 
observe that at a time of about 10 seconds after switch- 
off, when all but a few particles have merged with the 
cluster (Fig. [2]), the system seems to revert to a Haff-like 
behaviour. 

We see that there are no qualitative differences in the 
plots for different excitation methods. This implies that 
both methods produce a similar state of initial excitation. 



However, the experiments differ in sample density and 
particle speeds, which allows us to probe the parameter 
space and test the scaling of Haff's law with respect to 
these parameters. 

In conclusion, we have demonstrated a novel method 
to study granular cooling directly, validating Haff's law. 
The data can be described by the decreasing density of 
the free particles outside the cluster. This gives a descrip- 
tion of the cooling ground state of a granular gas even 
in the presence of clustering. Theoretically, the challenge 
remains to calculate the temporal dependence of the den- 
sity of free particles from clustering properties in order 
to have a full description of the system. 
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